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EVALUATION OF AN AEROELASTIC MODEL TECHNIQUE FOR 
PREDICTING AIRPLANE BUFFET LOADS 

By Perry W. Hanson 
Langley Research Center 

SUMMARY 

A wind-tunnel technique which makes use of a dynamically scaled aeroelastic model 
to predict full-scale airplane buffet loads during buffet boundary penetration is evaluated. 
A 1/8-scale flutter model of a fighter airplane with remotely controllable variable-sweep 
wings and trimming surfaces was used for the evaluation. The model was flown on. a 
cable-mount system which permitted high lift forces comparable to those in maneuvering 
flight. Bending moments and accelerations due to buffet were measured on the flutter 
model and compared with those measured on the full-scale airplane in an independent 
flight buffet research study. It is concluded that the technique can provide valuable infor- 
mation on airplane buffet load characteristics not available from any other source except 
flight test. 


INTRODUCTION 

Buffeting is the term loosely applied to the randomly fluctuating aerodynamic forces 
produced by separated flow or to the aeroelastic response of a structure to these forces. 
Buffeting of transport-type airplanes is of concern primarily with regard to onset condi- 
tions and shallow penetration of the compressibility buffet boundary due to upset condi- 
tions since this type airplane is not designed for high load factors. The concern is more 
for unpleasant vibrations and degradation of performance than for structural failure. 
Highly maneuverable high performance fighters, on the other hand, are subject to buffet 
flow due to shock — boundary-layer interaction, high angle of attack, and wake impinge- 
ment. They are typically flown well into the buffet boundary where the buffet loads and 
accelerations build up rapidly with penetration. It has generally been found that even for 
these conditions, wing buffet loads are not the critical consideration but rather (1) vibra- 
tions which subject fire control, navigation and reconnaissance equipment, instruments, 
and crew to a more severe operational environment and increase fatigue; (2) degradation 
of performance through increased drag and decreased lateral stability which detracts 
from tracking capability; and (3) excessive structural loads on tail control surfaces. With 
regard to item (3) the horizontal tail is normally considered to be the critical tail compo- 



nent for buffet loads. However, during recent wind-tunnel buffet studies on a fighter air- 
plane, it was found that the critical components at high angles of attack were the vertical 
tails which were vibrating primarily in a torsion mode. 

Although wing buffet phenomena have received much attention from investigators 
(refs. 1 to 8, for example), the ability to predict analytically the intensity of the buffet 
response as the buffet boundary is penetrated leaves much to be desired, and the ability 
to predict empennage buffet loads in turbulent wake flow is nonexistent. Although 
attempts are made during design to estimate buffet loads from past experience on sim- 
ilar configurations and from empirical means, the fact that buffet response phenomena 
are highly dependent on configuration detail and involve complex patterns of flow separa- 
tion has led to the use of wind-tunnel data as a primary tool in predicting buffet onset 
flight conditions and, to a lesser extent, buffet intensity levels. 

There are two fundamentally different methods of predicting the buffet intensity 
from wind-tunnel measurements on models of new configurations: (1) measure buffet 
accelerations or stresses on a dynamically scaled aeroelastic model and (2) measure the 
pressure fluctuations on a nominally rigid model and then calculate the dynamic response 
when these pressures act on the flexible structure. Each method has its advantages and 
disadvantages, but in either case there is a dearth of correlation information with which 
to evaluate the accuracy of the methods to predict full-scale buffet loads quantitatively, 
particularly for tail surfaces and transonic flight conditions. Most correlation studies 
have dealt chiefly with buffet onset predictions for wing buffet or with the relative magni- 
tude of buffet loads as the buffet boundary is penetrated (refs. 9 to 11, for example). Gen- 
erally, the predictions of buffet onset boundaries from wind-tunnel tests have been reason- 
ably successful (although the definitions of onset are quite varied). Load predictions from 
simplified models, however, have not been notably successful, with the possible fortuitous 
exception of wing buffet, and the literature reflects contradictory conclusions with regard 
to the several techniques. 

From considerations of the many factors involved, it would seem that the use of a 
complete dynamically and elastically scaled wind-tunnel model having a minimum of 
mount restraint would offer the potential for the most accurate prediction of absolute 
loads on both wing and tail surfaces and dynamic accelerations as the buffet boundary is 
penetrated. Models designed to provide freedom-from-flutter validation of new airplane 
designs closely satisfy the requirements for this type of buffet model. However, in prac- 
tice even with this technique there are factors that could possibly affect the accuracy of 
the data, such as Reynolds number, wind-tunnel-wall, and pitch-rate effects. 

An investigation was undertaken to develop and evaluate this wind-tunnel technique 
for predicting loads and accelerations due to buffet. The study included an assessment 
of (1) modeling requirements, (2) suitable means of approximating high load factor flight 
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conditions, and (3) the complexity of complete vehicle structural response to buffet flow. 
The method was evaluated by comparing the scaled buffet bending moments and accelera- 
tions measured on a 1/8-scale flutter model of the variable-sweep fighter airplane with 
those measured in a flight buffet research program (ref. 12). The investigation was con- 
ducted in the Langley transonic dynamics tunnel in the high subsonic and transonic speed 
ranges. 


SYMBOLS 

Physical quantities in this paper are given in both the International System of Units 
(SI) and U.S. Customary Units. The measurements and calculations were made in U.S. 
Customary Units. Factors relating the two systems are presented in reference 13. 

a root-mean-square normal acceleration, g units (lg = 9.80665 m/sec2) 

b characteristic length, m (ft) 

C/C cr ratio of generalized damping coefficient to critical value of damping 
Cj«qr normal-force coefficient, N/qS 

ACp rms root-mean-square fluctuating pressure coefficient, Ap/q 

root-mean-square fluctuation of pressure about the mean, N/m2 (lb/in^) 
bending stiffness, N-m2 (lb-in2) 

power spectrum of unsteady pressure coefficient, 1/Hz 
frequency, Hz 

nondimensional form of value of power spectrum of unsteady pressures at 
particular frequency 

torsional stiffness, N-m2 (Ib-in^) 

mass moment of inertia, m-N-sec2 (in-lb-sec2) 

aerodynamic scale factor, p r V r b r m r cu n>r 


Ap 

El 

F(f) 

f 

fm) 

GJ 

I 
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K a 

K ct 

k 
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m 

N 

q 

s 

V 

V 0 

A 

P 

a 

a> 

Subscripts: 


buffet acceleration scale factor, b r ^k r V2q rmr ^C'p/C cr j r ^ 

buffet bending-moment scale factor, b r ^k r ^ // ^q r ^Crpy , C C r) r ^ 
reduced frequency, bw/V 
Mach number 
mass, kg (slugs) 

total aerodynamic normal lifting force, N (lb) 
dynamic pressure, — pV 2 , N/m2 (lb/ft 2 ) 

Li 

wing panel reference area, 0.762 m2 (8.2 ft 2 ) for model 

free-stream velocity, m/sec (ft/sec) 

any particular nondimensional station 

sweep angle of leading edge of wing, deg 

airstream density, kg/m 2 (slugs/ft 2 ) 

root-mean-square bending moment, m-N (in-lb) 

circular frequency, 27rf , rad/sec 


A airplane 

a aerodynamic 

M model 

n nth natural vibration mode, n = 1, 2, 3, . . . 

r ratio, airplane to model 
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s structural 

T total 

t extraneous or "tare" value 

APPARATUS AND PROCEDURE 
Wind Tunnel 

This investigation was conducted in the Langley transonic dynamics tunnel. This 
tunnel has a 4.88-m-square (16-ft) test section with cropped corners and is a return-flow, 
variable-pressure, slotted-throat wind tunnel. It is capable of operation at stagnation 
pressures from about 1724 N/m2 (1/4 lb/in2) to slightly above atmospheric pressure and 
at Mach numbers up to 1.2. Mach number and dynamic pressure can be varied indepen- 
dently with either air or Freon- 12 used as a test medium. Freon- 12 was used in the 
present investigation. Stagnation temperature during the tests was approximately 48.8° C 
(120° F). 

An important factor in model dynamic measurements in wind tunnels is the airflow 
unsteadiness which, if excessive, may unduly affect the accuracy of the results or even 
mask the model dynamic characteristics to be evaluated. The ideal level of flow unstead- 
iness is, of course, zero. However, all tunnels have some flow unsteadiness and the 
acceptable level is determined by the type of dynamic measurements and the desirability 
of other tunnel operating characteristics. A measure of the suitability of this wind tunnel 
for buffet studies is depicted in figure 1. Figure 1(a) shows the variation of empty-tunnel 
sidewall rms fluctuating pressure coefficient with Mach number for the maximum dynamic 
pressure capability of the tunnel. The levels are of the order of fluctuating pressures in 
a turbulent flat-plate boundary layer (0.005 to 0.008 for the Mach number range considered 
(high subsonic and transonic)). (See ref. 14, for example.) Reference 15 presents a sum- 
mary of transonic wind-tunnel noise investigations that show the fluctuating pressures on 
the tunnel center line generally to be less than those at the walls. The fluctuating pres- 
sure levels shown in figure 1(a) may therefore be regarded as the probable upper limits 
likely to be encountered in the vicinity of the tunnel center line where the model was flown. 
These levels are generally quite low compared with the levels measured in the other tun- 
nels summarized in reference 15. It is of interest however to consider the power spectra 
of these pressure fluctuations. Figure 1(b) presents criteria for examining the suit- 
ability of a wind tunnel for buffet research, which take into account the frequency content 
of the unsteady tunnel flow. The criteria are based on the measurement of fluctuating 
static pressure on the tunnel wall and the results from studies of several different models 
in several different wind-tunnel test sections having varying degrees of flow unsteadiness 
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(ref. 16). In the figure, the ordinate yfF(f) is a nondimensional form of the value of the 
power spectrum at important natural frequencies of the model to be tested, where F(f) 
is the power spectrum of unsteady pressure divided by the dynamic pressure squared and 
f is the frequency in Hz. Along the left-hand side of the figure are shown the degrees of 
adequacy of the tunnel (from ref. 16) for two categories of buffet: light buffet generally 
pertaining to low-aspect-ratio, highly swept, thin surfaces and heavy buffet pertaining to 
relatively high-aspect-ratio, unswept, thick surfaces and surfaces in wake flow. From 
the figure it may be seen that the tunnel would be rated from good to adequate for heavy 
buffet and adequate to marginal for light buffet in the most critical area of M = 0.7 to 0.9. 
The horizontal tail would be the component most susceptible to these "extraneous" forces 
on the model in this range. 


Model and Support System 

Model description .- The model used in this investigation was a 1/8-size dynami- 
cally scaled aeroelastic model of a variable-sweep-wing fighter airplane that was initially 
designed and fabricated by the airplane manufacturer for use in flutter studies in the 
Langley transonic dynamics tunnel on a cable-mount system. The required stiffness was 
provided by steel spars in the wing and fuselage while ballasted segmented sections made 
from balsa and fiber glass provided the required shape and mass. The horizontal and 
vertical tails were of a stressed skin type construction (etched aluminum skin on balsa 
cores). The horizontal -tail pitch deflection and wing sweep angle were remotely control- 
lable and the horizontal tail could be differentially operated for roll control. Flow- 
through nacelles were incorporated, but mass flow through them was not simulated. The 
model was used successfully in flutter studies from the prototype phase to the present 
configuration and has undergone several updating changes. The only alteration made to 
the model for the buffet studies was the addition of strain gages and accelerometers at 
appropriate locations (see section "Instrumentation and Data Reduction") and a fitting for 
attaching the lift-balancing cable (see section "Lift-balancing device"). This model was 
designed to simulate a preprototype version of the fighter airplane and differs in some 
respects from the flight buffet test airplane. Although the model weight, center of gravity, 
and inertias compare favorably with the test airplane, the horizontal-tail pivot location on 
the test airplane is 0.127 m (5 in.) (full scale) forward of the model pivot location (less 
than 3-percent-root-chord difference). The model weighed 43.3 kg (95.5 lb) and had a 
full wing span of about 2.29 m (7.5 ft) at the 26° sweep position. Figure 2 is a photograph 
of the model suspended on the transonic dynamics tunnel cable-mount system while under- 
going flutter studies. The external stores and mounts shown were removed for the buffet 
tests. 

Model design scale factors .- The model was designed for a Mach number ratio of 1, 
and for a mass-density ratio scaling parameter of 1.0 at a full-scale airplane weight of 


6 



32 300 kg (71 200 lb), and a reduced velocity parameter of 0.837 (or, conversely, a 
reduced stiffness parameter of 1.195). The scale factors used in the model design were 
determined by wind-tunnel size, test medium, and tunnel performance capabilities and by 
the full-scale flight boundary. 

The present model was designed originally as a flutter model to demonstrate that 
the airplane would have the required flutter safety margin. Such models may be designed 
to be tested at wind-tunnel conditions that are scaled from the airplane actual flight bound- 
ary which has been increased by the flutter safety margin, in which case, k n r = 1; or, 
they may be designed with the model stiffness reduced in proportion to the flutter safety 
margin requirements and then tested at wind-tunnel conditions that are scaled values of 
the airplane actual flight boundary. 

The latter philosophy can alleviate practical problems in model fabrication by allow- 
ing the model to be less stiff for the same mass and was followed for the design of the 
present flutter model. Detailed discussion of dynamic aeroelastic model scaling theory 
is presented in reference 17. The basic design and derived scale factors are presented 
in table I. 

Model vibration characteristics .- Some of the more important model natural vibra- 
tion frequencies scaled to full-size airplane values are compared with airplane values in 
table n. These full-scale airplane values were measured on the first airplane by the 
manufacturer. The frequencies are almost the same as those given in reference 12 for 
the test airplane flown in the flight buffet tests. Model data were measured with the 
model suspended on very soft springs. The agreement between model and airplane values 
is considered to be very good. Although no attempt was made to control the model struc- 
tural damping, the table shows that the model and airplane structural damping ratios were 
generally comparable. 

Some of the model natural vibration mode shapes and frequencies of interest for the 
various wing sweep angles are shown in model scale in figure 3. All the vibration char- 
acteristics are given for the model ballasted to simulate the fuselage-full wing-empty 
fuel condition corresponding to an airplane weight of 32 300 kg (71 200 lb). The flight 
buffet data were measured on the test airplane with wing tanks empty, but the total air- 
plane weight varied from about 33 600 to 26 300 kg (74 000 to 58 000 lb) due to usage of 
fuel from the fuselage tank. This weight variation would have the primary effect of 
changing the natural frequencies slightly. 

Mount system .- For these tests the cable -mount system consisted of two cable 
loops. (See fig. 2.) A forward flying cable was reeved over pulleys in the model in the 
horizontal plane and secured to each tunnel sidewall about 9.14 m (30 ft) upstream of the 
model. The rear flying cable passed over pulleys in the model in the vertical plane and 
extended downstream of the model where one end was attached to the tunnel ceiling and 
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the other end was attached to a variable tensioning device through a soft spring. In addi- 
tion, a snubber cable system was provided for emergency restraint. It consisted of three 
cables which were fastened to the" model near the center of gravity and which- extended, out 
through the tunnel walls and ceiling to a shock absorber system and a remotely operable 
actuator. The snubber cables, shown taut in the photograph, are normally slack during 
data taking. The mount-system characteristics are described in more detail in refer- 
ence 18. 

The flexible umbilical is used to carry instrumentation leads and horizontal-tail and 
wing-sweep actuator motor power leads from the model. 

Lift-balancing device .- Normally the cable-mount system allows the model to fly 
rather freely near the tunnel center line with the model trimmed to support its own weight. 
In order to fly the model under high load factor conditions, it was necessary to counteract 
the lift that was in excess of the model weight; this was accomplished by means of a lift- 
balancing cable which was attached to the model near the center of gravity (c.g.) and which 
extended downward through the tunnel floor to a load link and then to a lift -balancing 
device which is illustrated schematically in figure 4. The cable is attached to a blowby, 
low-friction piston in a cylinder connected to large pressure reservoir tanks and to a 
large exhaust buffer tank. The large volume pressure reservoirs provide a relatively 
constant low level spring rate which allows the model to respond in rigid and elastic body 
dynamic motions with a minimum of restraint. The static lift -balancing force exerted by 
the lightweight piston is controlled remotely by varying the air pressure acting on the 
piston. The load link measures the force being exerted on the model by the device which 
is capable of producing up to 6670 N (1500 lb) of force. 

Instrumentation and Data Reduction 

For the buffet studies the model was instrumented to measure dynamic bending 
moments and accelerations at locations nearly the same as those on the airplane as 
possible. Four-active-arm, resistance-wire, temperature-compensated, strain-gage 
bridges to measure the left and right horizontal-tail bending moments at the root chords 
were bonded to the horizontal-tail pivot shafts at the same location as on the airplane. 
Strain-gage bridges also were bonded to the right and left wing spars and oriented to mea- 
sure bending moments at the same location as on the airplane (23.83 cm (9.38 in.), model 
scale, outboard of the wing pivot center measured along the 26-percent-chord line). Mea- 
surements were made only on the right wing of the airplane. An accelerometer to mea- 
sure normal dynamic accelerations was located in the model at the nominal c.g. location. 
The normal accelerations were measured on the airplane at a point 15.88 to 19.05 cm 
(6.25 to 7.5 in.), model scale, forward of this location, although in reference 12 they are 
referred to as being nominally at the center of gravity. In addition, an accelerometer 
was located in the model to measure the dynamic normal accelerations at the pilot station. 
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Other measurements made during the study were model angle of attack, force on the lift- 
balancing cable, left and right horizontal-tail deflection angles, wing sweep angles, and 
wind-tunnel parameters. 

Figure 5 is a schematic diagram of the signal conditioning and data-recording pro- 
cess. The signals from the model bending-moment strain gages and accelerometers 
were amplified by 3 -kHz carrier amplifiers. The amplified signals were recorded on 
oscillographs for total load monitoring. They were also routed through 200-Hz low-pass 
filters to eliminate frequencies above the range of interest and through dc blocking capac- 
itors to eliminate static bias to a 14-channel FM tape recorder. The dynamic portion of 
the signals was recorded at 9.53 cm/sec ^3| in./secj for later electronic analysis. The 
analog signals were converted to engineering units and digitized at a sampling rate of 
500 samples/sec. Standard digital -computer random-process analysis techniques were 
used to generate power spectral densities and rms levels of the model response. The 
model data were analyzed over a frequency range from 3 to 150 Hz, with a frequency band 
width of 3 Hz. (The airplane data of ref. 12 were analyzed to 50 Hz with a frequency band 
width of 1 Hz.) The model data sample periods were approximately 30 sec, whereas the 
airplane data sample periods were 3 sec. Signals from the horizontal-tail and wing posi- 
tion indicators, model angle of attack servo accelerometer, lift-balancing cable load link, 
and tunnel pressure and temperature sensors were routed through digitizing potentiom- 
eters to a card punch. Test-point identification information served to correlate the tape- 
recorder data with the punch-card data. 

From repeatability tests it is estimated that the maximum error of the wind-on 
measurements of model normal-force coefficients, dynamic bending moments, and 
dynamic accelerations was approximately ±5 percent; maximum error in angle of attack 
was approximately ±1/2°. Tunnel flow parameters are generally accurate to less than 
±1 percent except that the maximum error in Mach number is about ±1/2 percent. 

Test Conditions and Procedures 

Measurements were made on the model at tunnel conditions (Mach number and 
dynamic pressure) that simulated most of the full-scale flight conditions during the flight 
buffet measurements reported in reference 12. In the wind-tunnel tests the procedure 
was to stabilize the tunnel at a Mach number and dynamic pressure scaled from averaged 
values measured during the particular flight buffet test maneuver. The model was first 
flown at a given wing sweep angle on the cable -mount system just lifting its own weight. 
The horizontal tails were then positioned to produce discrete increments in model angle 
of attack. At each selected angle of attack a 30-sec tape data sample of model bending 
moments and accelerations was recorded, along with the tail positions, angle of attack, 
lift-balancing cable load, and tunnel conditions. The angle of attack was increased in 
increments until no further increase was possible either because maximum tail deflection 
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angles or maximum model design loads were reached, or because of model instability. 
The wing sweep angle was changed to a new setting and the procedure was repeated for 
the other sweep angles. 

The model and airplane relative buffet onset boundaries and buffet boundary pene- 
trations are compared in figure 6. The reasons for termination of the model penetra- 
tions are indicated in the figure. 

Figure 7 compares the variations of the model normal-force coefficients with angle 
of attack with those of the airplane for various Mach numbers and wing sweep angles. 

For comparable test conditions the flight Reynolds numbers (based on the wing mean 
chord at 16° sweep) varied from 2.2 x 10® to 3.1 x 10®, whereas the model Reynolds num- 
bers varied from 0.77 x 10® to 1.18 x 10®. The model and airplane data are felt to com- 
pare favorably. 


DATA ANALYSIS 

Scaling Buffet Loads and Accelerations 

For wind-tunnel buffet studies on models to be useful in predicting full-scale buffet 
loads, the necessary buffet response scaling relationships must be known. Over the years 
investigators have developed buffet scaling relationships that generally have been much 
simplified and oriented to the particular study at hand (refs. 4 and 5, for example). The 
complexity of these relationships depends on the complexity of the dynamic system 
involved. For a complete airplane structure the system is indeed quite complex and 
certain simplifying assumptions are required, such as considering components of the 
structure individually (wing, horizontal tail, etc.) in order to derive workable relation- 
ships. The development of buffet loads scaling relationships normally follows the method 
of generalized harmonic analysis. The method is treated in some length in reference 19 
and was first applied to the analysis of buffeting in reference 1. One of the more general 
developments which considers distributed load and structural characteristics and dis- 
crete modal contributions to the total response is presented in reference 20. The devel- 
opment requires that the model frequency spectrum be reasonably representative of the 
full-scale vehicle, or at least be representative of the important vibration modes involved 
in the vehicle buffet response. Although the scaling relationships developed in refer- 
ence 20 were derived for scaling launch vehicle buffet responses, they are applicable to 
any linear multi -degree -of-freedom system. They are used as a basis for buffet loads 
scaling in this study. 

From equation (15) of the appendix of reference 20 the buffet bending-moment scal- 
ing relationship may be written in the present notation as 
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WJ 

1 k ", 


5Wn,M \ C cr/n,Mr^ 


C cr/ n ,A \ C cr/n,A_ 


[CL,nW]/ 2 K) n ,, 


o(n Q ) T ^ airplane total root-mean-square buffet bending moment at a particular location 
b r airplane-to-model length ratio 

q r airplane-to-model dynamic pressure ratio 

k n r airplane-to-model reduced frequency ratio for nth natural vibration mode, 

br w n,r/v r 

Cl n (^n) airplane-to-model ratio of power spectrum of effective random aero- 

r dynamic lift coefficient 

a^ 7 °)n model root-mean-square buffet bending moment in nth vibration mode at a 
particular location 


Ca_\ + (Cs 


sum of model aerodynamic and structural damping in nth vibra- 
tion mode 


Ccr /n,A V C cr/ n;A 


sum of airplane aerodynamic and structural damping in nth 
vibration mode 


An expression similar to equation (1) can be derived for the airplane root-mean - 
square buffet acceleration at a particular location 
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where 


m r airplane -to -model mass ratio - 

a (m.) model root-mean-square buffet acceleration in nth vibration mode at partic- 

' u 'n,M 

ular location 

From equations (1) and (2) it can be seen that the total mean-square bending moment 
and acceleration is a superposition of single-degree-of-freedom responses, that is, the 
derivation treats each mode as a separate independent system with modal coupling 
neglected. Thus, the model/airplane scaling relationships are dependent on the struc- 
tural and aerodynamic damping ratio, the reduced-frequency ratio of each vibration mode 

of interest, and the term n(^n)J ■ Ideally, the reduced-frequency ratio should be 

unity in order that the effective aerodynamic lift coefficients of the model and airplane be 
identical. 

Because of slight variations of the speed of sound in the wind tunnel with respect to 
the speed of sound in flight (the airplane tests were conducted at several altitudes) and 
because the present model was designed to demonstrate the required flutter margin by a 
reduction in stiffness, the actual reduced-frequency ratios varied from the optimum value 
of one by 25 to 37 percent. The actual values of the square root of the reduced-frequency 
ratios are shown in table III. These values are based on a frequency ratio u> r approx- 
imately equal to 1/3 for all modes (determined from measured airplane and model fre- 
quencies). For the range of reduced frequencies of this study, these variations of k r 

are considered to be negligible and the value of the term n (kj^j is assumed to be 
unity. _____ - 1 

For dynamically scaled aeroelastic models with relatively low structural damping, 
the damping term in brackets in equations (1) and (2) will normally be very near unity. 

(See section "Modeling Considerations.") However, it was of interest to estimate the 
value of the term in the present study. The model structural damping (c s /C cr ) n and 

airplane structural damping (c S/ /C cr ) n ^ were measured for the modes considered to 

be of primary importance in the buffet response. The model aerodynamic damping 
(Ca/ C cr) n ^ was calculated based on reference 21 since the flutter model used did not 

have sufficient instrumentation to measure aerodynamic damping in buffet flow. Calcu- 
lated values of airplane aerodynamic damping were also used since no measured values 
were available from the flight buffet tests. The airplane aerodynamic damping can be 
related to the model aerodynamic damping by 
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( 3 ) 


Ca\ _ / Pr^r^r^ j / ^a\ 
£ ct L,A \ mrWn ’ r f\ Ccr )n,M 


(see ref. 21, for example) or 



(4) 


where 


K D = 


P r V r b r 

m r a) n,r 


During flight flutter tests at A = 26°, some limited measurements of the airplane aero- 
dynamic damping in the wing bending mode were made. The measured airplane values of 
0.77 to 0.82 agreed well with the calculated values of 0.75 to 0.94. 

Since the frequency ratio is taken to be the same value for all modes of interest 
(u ) T •= 1/3^, the reduced-frequency ratio will be the same for all modes and equations (1) 
and (2) may be written as 


ct 2 Wt,A = b r 6( lr 2k r 



(5) 


and 


OO 

a 2 f?? 1 ) = b r ^q r 2 k r — — Y — a 2 (r/ ) 

[ °>T,A r mr 2 L x (G T /C cr ) r n V<«, 


where 


(6) 
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Actually it is not generally possible to measure directly the bending -moment 
a (^°)n M ° r acce * era ^ on a (^o) n buffet response of the model in a particular mode. 

Rather, the total response at a particular location is measured and the "modal com- 

position" of the indicated response is estimated from power spectra. Several conditions 
are inferred in this application of the scaling relationships: (1) The "natural vibration 
modes" of importance in the total response can be identified, are well separated, and are 
lowly damped; (2) the total measured model response can be treated as a summation of 


OO 


individual modal responses, that is, 



2 


n,M : 


and (3) the structural and aerody- 


n=l 

namic damping of the modes are known or can be estimated. Although condition (1) may 
not be strictly met for the horizontal tail and c.g. acceleration responses because in some 
cases the modes were not well separated, it was assumed that these responses came 
acceptably close to meeting this condition for application of the scaling relationships. 
Therefore, the airplane buffet response at a particular location can be determined by the 
modal summations indicated in equations (5) and (6). In practice, this may prove to be a 
difficult process. In the present case, the character of the responses and the fact that 
the model was a good dynamic representation of the airplane permitted some simplifica- 
tions. The measured total wing bending-moment response (of both the model and the air- 
plane) was very similar to that of a single-degree-of-freedom system with the response 
frequency equal to the wing first bending frequency. (See fig. 8(a).) In contrast, power 
spectra of the horizontal tail and c.g. accelerations (figs. 8(b) and 8(c)) reveal some 
energy at frequencies approximating those of several natural vibration modes. However, 
the model spectra generally indicated that the major portion of the total response was 
near frequencies corresponding to the horizontal tail bending, fuselage vertical bending, 
and wing first bending modes. Since the total damping ratios (cj^C cr j^ ^ for each of 


these modes was very nearly the same (table III), it was expedient to assume that the 
damping associated with all the modes involved in the horizontal-tail total bending-moment 
response was that of the horizontal -tail bending mode. Similarly, the total damping ratio 
associated with the total c.g. acceleration response was considered to be that of the fuse- 
lage vertical bending mode. 


Equations (5) and (6) are then simply 


CT ( r? °) T ,A - a A ~ K ct ct M 

and 

a W T)A - a A = K a a M 


(7) 


( 8 ) 
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where 


K 


CT 


= b r ^k r 



K a = b r 


2. 1/2 


q r m r 



-1 / C T 


(9) 


(10) 


and ct m is the total rms buffet bending moment measured on the model wing or horizon- 
tal tail, ajyj is the total rms buffet acceleration measured on the model at the c.g. and 
pilot station, and a ^ and a^ are the corresponding full-scale bending moments and 
accelerations. The total damping ratio ^C>j>^/C cr j is that of the wing first bending 

mode, the horizontal-tail bending mode, or the fuselage vertical bending mode for wing, 
horizontal tail, or c.g. acceleration measurements, respectively. The values of the flow- 
related scaling ratios and scaling constants used in converting from model to full-scale 
values are shown in table III. 

Extraction of Extraneous Response 

Both the flight- and wind-tunnel-model dynamic loads and accelerations, as mea- 
sured, are the result of both buffet and extraneous forces. The airplane instrumentation 
measures structural response to atmospheric turbulence, sudden actuator movements, 
and engine and auxiliary motor vibrations, in addition to the response of the structure to 
buffet flow. The model instrumentation measures response to tunnel turbulence and 
mounting cable vibrations in addition to the response to buffet flow. Generally these 
extraneous responses are relatively constant for a particular Mach number and dynamic 
pressure as angle of attack or normal force is varied and can be considered "tare" values 
to be subtracted from the total response. Figure 9 illustrates the manner in which the 
model tare values were determined. The total measured model and airplane dynamic 
bending moments and accelerations are shown as a function of normal-force coefficient. 
For comparison purposes, the total model response has been scaled to full-scale values. 
Figure 9(a) is typical of the model response for the 26° wing sweep condition for all Mach 
numbers and for the higher sweep angles at low Mach numbers. Figure 9(b) is more gen- 
erally representative of the response at higher wing sweeps at the higher Mach numbers. 
The wing bending moments and c.g. accelerations are seen to increase only slightly until 
a certain value of is reached, after which the level of response increases rapidly. 

This characteristic also is true for the horizontal-tail bending moments for the low wing 
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sweep conditions. Model buffet response is said to begin at the intersection of straight 
lines faired through the total response prior to and after buffet onset. The values of the 
accelerations or bending moments-at this- point are taken to bethe tare values a^ or 
ctj., respectively. For those instances where the horizontal -tail response did not exhibit 
a clear change in character, the tare values were determined by multiplying the relatively 
well-defined tare values at M = 0.52 by the ratio of the tunnel turbulence level at the 
Mach number in question to the level at M = 0.52 (as determined from fig. 1(b)). In 
most instances the tare values determined in this manner were approximately the same 
as the minimum values of bending -moment response observed over the range of Cjg 
variation. 


If there is no correlation between tunnel unsteadiness and the buffet response the 
actual model root-mean-square buffet bending moments are given by 


a 


M “ 




( 11 ) 


where Op is the measured total rms bending moments. 

Similarly, the model buffet accelerations a^ are given by 


a M " 




( 12 ) 


where a-p is the measured total rms accelerations. 

The measured airplane rms bending moments and accelerations given in refer- 
ence 12 and shown normalized by airplane design values in figure 9 are total values. For 
the flight buffet studies extraneous nonbuffet related structural responses were considered 
to be insignificant since the flight data were purposely taken under minimal turbulence con- 
ditions. From unpublished data supplied by the authors of reference 12 and considered by 
them to be typical of prebuffet levels for all the flight data, it was concluded that for con- 
sistency in the model/flight comparisons the following small tare values would be 
appropriate: 


Measurement 

Prebuffet tare values, 
percent design 

rms wing bending moment (right) 

0.112 

rms horizontal -tail bending moment: 


Left 

0.147 

Right 

0.158 

rms acceleration at c.g 

0.177 
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These nonbuffet related values were subtracted from the total airplane values presented 
in reference 12 in the same manner as for the model data. Buffet onset on the airplane 
was said to occur when the peak-to-peak c.g. accelerations exceeded approximately 
0.8 percent of the design load factor. 

It may be noted that the mathematical procedure required to separate the nonbuffet 
related response from the total (taking the square root of the difference of squares) 
results in tare values having a minimum effect at the higher values of total response 
where the greater interest lies for predicting buffet loads. 

Figure 9 also illustrates an unfortunate aspect of the flight buffet studies with 
regard to comparative evaluations, namely the low level of the buffet loads and accelera- 
tions due to the previously mentioned load-factor and angle-of-attack limitations that had 
been imposed on the preprototype airplane. The maximum dynamic component of the 
loads measured was generally not greater than 1 percent of design loads and the maxi- 
mum acceleration response was generally little greater than 2 percent of design load fac- 
tor. These low levels result in a rather severe test of the buffet loads prediction cap- 
ability of the present method due to the relatively larger values of extraneous response 
and effects over which the investigator has little or no control. 

DISCUSSION OF AEROELASTIC MODEL TECHNIQUE 

Some of the more important modeling considerations, including compromises, nec- 
essary to evaluate effects of buffet flow are discussed and related to the buffet loads scal- 
ing relations. Model testing technique, including limitations, deficiencies, and potential 
impact on scaled results are then discussed. 

Modeling Considerations 

The goal, of course, is to subject a reduced size model that accurately reflects the 
shape and structural characteristics of the airplane to wind-tunnel flow that accurately 
simulates the airflow over the airplane and to support the model in the wind tunnel in a 
manner such that the model and airplane degrees of freedom and inertia forces are prop- 
erly related. In practice, it is impossible to meet all these requirements, but a complete 
dynamically scaled aeroelastic model suspended on the cable -mount system described 
earlier comes reasonably close to meeting them. The model should be as large as pos- 
sible, consistent with the tunnel test section size to facilitate model fabrication and, for 
buffet models, to provide a Reynolds number as large as possible. Buffet flow is influ- 
enced by the nature of the boundary layer which, of course, can be a function of Reynolds 
number. Because of wind-tunnel performance limitations, practical model fabrication 
considerations, and the overriding importance of testing the model at a mass-density 
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ratio comparable to that of the airplane, no attempt normally is made to simulate 
Reynolds number in high-speed flutter models. The same compromise would have to be 
made for high-speed buffet models. Another factor affecting the boundary layer which 
may be more important for buffet models than for flutter models is the skin surface 
smoothness. The model used in the present studies was of a "segmented pod” type con- 
struction which results in small surface discontinuities at the edges of the individual seg- 
ments. However, the Reynolds number scale effects and surface roughness are believed 
primarily to influence the buffet onset boundary and residual model response prior to buf- 
fet onset rather than the load intensities once the buffet boundary has been penetrated. 
Figure 6 indicates that even the buffet onset conditions were predicted reasonably well 
except for the 72° wing sweep condition and figure 7 shows that the variation of normal- 
force coefficient with angle of attack was essentially the same for the airplane and model. 
Although the locations of local shocks and commencement of separated flow may have 
been different for the model and airplane, their integrated effects on the structural 
response apparently were small. 

Once the model size is selected, the stiffness and mass scale ratios remain to be 
determined. These ratios are a function of model size. The airplane flight conditions 
must be simulated (Mach number and altitude) within the wind-tunnel operational capabil- 
ities. For a Mach number ratio of 1 and a nominal ratio of airplane to tunnel speed of 
sound, the resulting density and velocity ratios determine the nominal mass and stiffness 
ratios. The model is then fabricated to obtain the scaled mass and stiffness distributions 
as closely as possible. The natural vibration mode shapes of the successful model will be 
very nearly the same as the airplane. If the measured node lines of the vibration modes 
of the model and the calculated or measured node lines of the airplane are very near the 
same location and if the scaled model frequencies are within a few percent of the desired 
value, the model is said to be satisfactory. 

From equations (1) and (2) it is clear that the relationship between the total damping 
(structural plus aerodynamic) is as important in scaling buffet response as the stiffness 
and mass characteristics. For a properly scaled model tested at exactly the scaled flight 

conditions, the airplane aerodynamic damping (C a /C cr j is related to the model aero- 

' ' 'n,A 

dynamic damping (C a /C cr ) by the reciprocal of the reduced stiffness ratio k r ; 
v ' ' n,M 

therefore when k r = 1, they are equal. (The aerodynamic damping scale factor Kjq 
of eq. (4) differs from l/k r to account for off-design test conditions of the model 
used in the study.) If it could be arranged to have the model structural damping equal to 
that of the airplane for the various vibration modes, then the total damping ratio would be 
equal to 1. Unfortunately, modeling technology has not advanced to the point where this is 
possible and attempts to control model structural damping are usually limited to keeping 
the model structural damping to reasonably small fractions of critical damping or pro- 
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viding an artificial means of varying the structural damping over a limited range. This 
latter method was used in the damping studies of reference 22 to confirm the relationship 
between buffet bending moment and total damping coefficient 


(C T /C cr ) 1/2 

The measured structural damping coefficients of the model used in the present study were 
approximately the same as the airplane structural damping coefficients for the primary 
modes of interest (table II), and were roughly an order of magnitude smaller than the cal- 
culated model aerodynamic damping (table IE). Since the aerodynamic scaling factor Kj) 
did not vary by more than +5 percent and -14 percent from 1.0 for all the test conditions, 
the resulting total damping ratio for the present model has only a small effect on the bend- 
ing moment and acceleration scaling factors K CT and Ka. This situation is in contrast 
to the one where solid metal aerodynamic force models are used to estimate wing buffet 
loads. Although in some instances such models may have wing reduced frequencies 
approximately equal to those of the airplane when tested at suitable velocities, the model 
aerodynamic damping Ca/C cr will be low due to the relatively large model mass 
(c cr = 2mo)j so that the model structural damping coefficient is likely to predominate, 
whereas the airplane aerodynamic damping coefficient is likely to predominate except for 
flight conditions approaching flutter. 

Even with the use of dynamically scaled aeroelastic models, certain simplifying 
assumptions are usually made. For the present studies, it was assumed that wind-on 
mode shapes were essentially the same as wind-off shapes and that there were no signif- 
icant differences in model and airplane mode shapes that were important in the buffet 
response. It was also assumed that the model and airplane structural damping was con- 
stant; that is, the structural damping in a particular mode was independent of vibration 
amplitude, temperature, and flow conditions. In fact, the structural damping may vary 
with amplitude and temperature so that the accuracy of the method could be improved by 
incorporating in the model a means for measuring the total wind-on damping in the vibra- 
tion modes of interest or by determining the total damping of the model response by con- 
sideration of the half-power point bandwidth of resonant peaks in the power spectra. (See 
ref. 23, for example.) If damping is to be determined in this manner, appropriate sam- 
pling rates and filter bandwidths must be used and spectrum analyzer characteristics 
must be accounted for. The technique was not used in the present study because the 
power spectra parameters (filter bandwidth, frequency range, etc.) were chosen to be 
compatible with the flight spectra for comparison purposes rather than for accuracy in 
determining model damping. 
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Testing Technique 

A basic difference in the flight and wind-tunnel buffet studies was the manner in 
which the structure was introduced to buffeting flow. The model angle of attack was 
slowly increased at essentially constant flow conditions and the response measurements 
were made at a constant angle of attack, whereas the airplane approached buffeting in 
maneuvering flight with a relatively rapid change in angle of attack through buffet onset 
and penetration with some variation in Mach number and altitude. The limits for the 
flight test conditions were either a maximum angle of attack of 18° or a maximum load 
factor of approximately 3g. The tunnel flow conditions were set to simulate the airplane 
average Mach number and dynamic pressure during the higher load, postbuffet portions of 
the maneuvers. However, it was not possible to change the model angle of attack at the 
rates experienced by the airplane (generally 1 to 5 deg/sec). The horizontal tails on both 
the airplane and model are used for pitch and roll control. However, on the model they 
are used as trim rather than maneuver devices and, hence, have a very slow deflection 
rate. Several studies have indicated that the maximum lift force attainable by a wing in 
unsteady motion increases with its pitch rate, that pitch rate tends to delay buffet onset 
and to increase the magnitude of the buffet loads. For high pitch rates such as might be 
encountered in helicopter or compressor blades the effect can be substantial, but for pitch 
rates comparable to those attainable by airplanes, even fighter types, these studies indi- 
cate the effect is generally less than a 10-percent increase in maximum lift coefficient 
for the Mach number range covered by the present studies. (See refs. 4 and 24, for 
example.) Evidence also suggests that the maximum buffet load at a given penetration 
into the buffet boundary tends to increase with time to some asymptotic value. The 
model, of course, can be subjected to the buffet flow for considerably longer periods than 
the maneuvering airplane, so that the pitch rate and time effects may tend to cancel. 

The test technique used in the present buffet studies exhibited a deficiency in 
horizontal-tail load simulation that has since been minimized; that is, the deflection of 
the horizontal tails relative to the model angle of attack was not the same for the model 
and airplane. On the airplane, except for trimmed flight conditions, the instantaneous 
deflection angle of the horizontal tails relative to the instantaneous angle of attack for a 
particular dynamic pressure will vary, depending on the type of maneuver (rolling or 
straight pull-up), the suddenness of the tail deflection, and whether the deflection is to 
initiate or terminate the maneuver. On the model, however, the horizontal tails were 
used to statically trim the model to the required angle of attack. In addition, the model 
horizontal tails had to overcome a pitching moment that the full-scale tails did not have 
to contend with (the resisting moment of the cable -mount system) so that the effective 
tail power of the model was less than that of the airplane. A means has since been devel- 
oped that allows horizontal-tail deflections independent of model angle of attack within 
certain limitations. A thin cable attached to the model forward of the effective pitch axis 
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extends vertically to a pitch-moment device located overhead outside the test section that 
is quite similar to the lift-balancing device described earlier. A load link is incorporated 
in the cable to allow determination of the pitching moment and lift that the device exerts 
on the model. With this system the tail buffet loads can be measured as a function of 
horizontal-tail deflection angle over a wide range of values for various discrete model 
angles of attack. Any effect of the cable on the model structural response can be mini- 
mized by attaching the cable near the fuselage vertical bending mode node line. 

DISCUSSION OF RESULTS 
Characteristics of Buffet Response 

Power spectral density plots of the model structural response (bending moments 
and accelerations) were used to analyze the character of the buffet loads for the various 
flight conditions studied. Almost 500 plots were generated from which some conclusions 
can be made regarding the modal contributions to the total response. Many were repeti- 
tious and only samples or examples which serve to illustrate a point are presented in this 
report. Since the comparisons of model and airplane buffet spectra indicate that the 
model buffet characteristics were reasonably similar to those of the airplane, the power 
spectra of the various measurements made on the model were studied to determine the 
effect of wing sweep and Mach number on the modal contributions to the total response. 
Although there were some inconsistencies or anomalies, particularly in the horizontal- 
tail data, some general observations can be made, bearing in mind that because of the 
close proximity of frequencies of several of the natural vibration modes and because 

these frequencies may change some with flow velocity, it is sometimes difficult to ascribe 

1 

a particular natural vibration mode to a particular "spike" on the power spectral density 
plot. In general, the modal composition of the total buffet response was dependent on the 
wing sweep /angle, Mach number, and in some cases, depth of penetration into the buffet 
region. 

Figure 8 shows some sample comparisons of model and airplane response spectra 
considered to be typical. The data are for a normal-force coefficient near that for max- 
imum buffet. 

Wing bending moment . - In figure 8(a) the right-wing bending moments of the model 
and airplane, normalized by the root-mean-square bending moments, are shown in terms 
of full-scale frequencies. The airplane data are from reference 12. The model and air- 
plane wing bending-moment responses are seen to be almost the same and apparently 
consist solely of the wing first symmetrical bending mode. Measurements of wing nor- 
mal shear force on the airplane indicated a slight contribution of the wing second bending 
mode to the total wing shear response. The model was not instrumented to measure this 
quantity. The peak response occurred at a frequency that correlated well with the wing 


21 



first symmetrical bending mode for all sweep angles and Mach numbers tested, both prior 
to buffet onset and during buffet. At very high values of normal-force coefficients, very 
slight amounts of energy were indicated at frequencies corresponding roughly to wing first 
torsion and/or horizontal-tail pitch. 

Horizontal-tail bending moment .- Figure 8(b) compares the model horizontal-tail 
root bending-moment spectrum with unpublished airplane data for the left horizontal tail. 
Here the model response is primarily in the horizontal -tail first symmetrical bending 
mode with secondary response in the fuselage vertical bending mode. Unpublished air- 
plane spectra indicate that the primary response varied between horizontal-tail bending 
and fuselage vertical bending depending on flight condition and whether the right or left' 
tail was being considered. Some of the variance between airplane and model response 
may be due to the lack of simulation in some areas mentioned previously (slightly different 
pivot axis locations for the model and airplane and different tail deflection angles relative 
to the model and airplane angles of attack). At other Mach numbers and wing sweep angles 
the primary horizontal-tail bending -moment response also occurred at frequencies corre- 
sponding to the horizontal -tail symmetrical first bending frequency. Secondary responses 
(spectra peaks less than half the primary response peaks) apparently occurred in the 
horizontal-tail pitch, wing first torsion, wing first bending, and fuselage vertical bending 
modes with relative magnitudes dependent on Mach number, wing sweep angle, angle of 
attack, and horizontal-tail deflection angle. Since the tail deflection angle, angle of attack, 
and wing sweep angle could not be varied independently for these studies (see section 
"Testing Technique"), it was not possible to separate the effects of these parameters. 

Center-of-gravity acceleration .- Figure 8(c) compares the acceleration spectrum of 
the model at the center of gravity with an unpublished airplane spectrum measured near 
the airplane nominal center of gravity for corresponding flow conditions. In this compari- 
son the airplane spectrum reveals a large amount of energy, around 40 hertz, not indicated 
in the model spectrum. The source of this response is not known. It may be a local res- 
onance of the structural member on which the accelerometer was mounted or, less likely, 
from the horizontal-tail antisymmetric torsion mode. Other available unpublished air- 
plane acceleration spectra indicate that at other sweep angles this energy peak is much 
suppressed or not apparent. In these instances the primary airplane c.g. acceleration 
response was in the fuselage vertical bending mode as was the model response. The pri- 
mary model response at all other wing sweeps and test Mach numbers was in the fuselage 
vertical bending mode. Secondary responses occurred in the horizontal-tail pitch mode 
with a limited amount of horizontal-tail fore-and-aft bending and wing first torsion modes. 
The relative magnitudes of the secondary responses were not as dependent on Mach num- 
ber as was the case for pilot station accelerations. 
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Pilot station accelerations .- The model buffet studies included measurement of the 
acceleration response at the pilot station, as well as at the center of gravity. Figure 10, 
which compares the acceleration spectra of measurements made at the pilot station and 
center of gravity at the same time, illustrates that the character of buffet accelerations 
is dependent on the measurement location. At the pilot station the buffet acceleration 
response primarily consists of the wing symmetric second bending mode, whereas at the 
center of gravity, the primary response is in the fuselage vertical bending mode with a 
significant amount of energy in the horizontal-tail pitch mode. Such differences in the 
frequency content of the response are particularly significant in considerations of the 
acceleration environment to which sensitive avionics equipment may be subjected. It 
should be noted also that since the acceleration is proportional to the frequency squared, 
high frequency structural response can cause quite high rms accelerations even though 
the amplitude or stress level is low. The fallacy of specifying some limit magnitude of 
buffet acceleration at an arbitrary location, such as the center of gravity, without consid- 
ering the primary frequencies involved in the response, their relationship to location on 
the structure, and the important resonant vibration characteristics of the devices is 
obvious. 

Prior to buffet onset, or at low levels of buffet, energy appeared at several frequen- 
cies about equally, but as the buffet condition became more fully developed, one mode gen- 
erally would become predominant. The primary response at all wing sweeps and Mach 
numbers was at a frequency corresponding to the wing second symmetrical bending mode. 
Secondary responses occurred at wing first bending, fuselage vertical bending, and some 
at the horizontal-tail fore-and-aft bending frequencies with the relative magnitudes 
depending on Mach number and wing sweep angle. At 26° sweep, the major secondary 
response was at the wing first symmetrical bending frequency with the amount of wing 
first bending energy relative to the primary response in wing second bending increasing 
with Mach number. At 50° and 72° sweep, the secondary response was at the fuselage 
vertical bending frequency instead of the wing first bending frequency as was the case for 
a wing sweep of 26°. 

Figure 11 compares the magnitudes of the buffet rms accelerations measured on the 
model at the center of gravity and at the pilot station. The figure shows that the manner 
in which the buffet accelerations change with normal-force coefficient generally is the 
same at both locations and is characterized by a sharp increase in magnitude with 
increase in normal-force coefficient for all Mach numbers and wing sweep angles except 
at M = 0.8 and M = 1.17 at A = 72°. For these conditions, the buffet accelerations 
at both locations remained at a relatively high constant level for all values of normal- 
force coefficient, a characteristic exhibited by the airplane at the higher Mach number at 
72° sweep. The intensity of the buffet accelerations at the pilot station are seen to be 
generally less than at the center of gravity. 
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Comparisons of Scaled Model and Airplane Buffet 
Bending Moments and Accelerations 

Power spectral density comparisons have shown' that the model and airplane buffet 
response characteristics were reasonably similar. The root-mean-square levels of the 
airplane and buffet bending moments and accelerations predicted from model data normal- 
ized by airplane design values, are compared in figure 12. The degree of correlation 
ranges from good to bad, depending on wing sweep angle, Mach number, and type of mea- 
surement. Except as noted, buffet bending moments are shown for both the right and left 
wings on the model. Only the right airplane wing was instrumented. It is interesting to 
note that although buffet onset occurred on both model wings at the same value of normal- 
force coefficient, the right wing consistently exhibited a larger response than the left. In 
addition, the tare response of the right wing was consistently higher than for the left wing. 
(See fig. 9.) A possible explanation is that the right wing experienced an effective higher 
angle of attack than the left wing, either due to differential twist in the wings or due to dif- 
ferential flow angularity. During tests the model tended to roll to the left so that a less 
negative horizontal-tail deflection angle was consistently required for the left tail to keep 
the model in level flight. 

Before assessing the comparisons shown in figure 12 it should be mentioned that for 
the flight data the Mach number decreases with increasing normal-force coefficient with 
most of the Mach number change occurring between the two largest values of normal-force 
coefficient. This decreasing Mach number during the maneuver may have an alleviating 
effect on the flight buffet response at the higher normal-force coefficients. 

A = 26° .- Figure 12(a) compares the wing and horizontal-tail buffet bending moments 
and the c,g. buffet accelerations for a wing sweep of 26° and a nominal Mach number of 
0.52. The airplane Mach number varied from 0.53 at the lowest buffet normal-force coef- 
ficient to 0.49 at the highest whereas the model data were for M = 0.52. After buffet onset, 
the bending moments on both model wings increase rapidly with normal-force coefficient 
and compare favorably with the flight data. The apparent decrease in flight buffet bending 
moment indicated by the data point at the maximum may be attributable to the 

decreasing Mach number or to normal scatter in the data. This amount of scatter would 
not be unusual for such flight measurements (ref. 4, for example) since the necessarily 
short data samples may not be long enough to exhibit completely stationary characteristics. 

The airplane and model horizontal-tail buffet bending moments are seen to compare 

I 

very well for the test condition of A = 26°, M = 0.52. There was no consistency in the 
maximum bending moments occurring on a particular side, either on the model or the 
airplane. 
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The airplane c.g. buffet acceleration does not show the rapid increase in magnitude 
with normal-force coefficient exhibited by the model data. At M = 0.76 (fig. 12(b)), the 
decrease in airplane Mach number with increasing is minimal (0.76 to 0.75) and for 

this case the model and airplane c.g. buffet accelerations both exhibit a rapid increase in 
intensity with increasing Cjj and agree very well. The airplane and model wing and tail 
bending moments also are in good agreement. 

Figure 12(c) compares model and airplane data for a nominal Mach number of 0.81 
at a wing sweep of 26°. The airplane Mach number varied from 0.83 to 0.81. All three 
measurements are seen to correlate very well. 

A = 50° .- At M = 0.52 and A = 50° (fig. 12(d)), the airplane and model wing buf- 
fet bending moment are reasonably comparable over the Cjg range duplicated. Values 
of Cn greater than those indicated for the model were not attainable due to insufficient 
tail power as discussed earlier. The model right horizontal -tail bending moments are 
seen to increase much more rapidly than the model left tail moments or the airplane bend- 
ing moments. The reader is reminded that the model horizontal -tail deflections normally 
were not the same as those of the airplane for equivalent normal-force coefficients. The 
model c.g. buffet accelerations are seen to increase almost linearly with whereas 
the rate of increase for the airplane apparently decreases at the higher normal -force 
coefficients; this is also true at M M = 0.80 (fig. 12(e)). The wing buffet bending 
moments predicted from model data tend to be slightly high as compared with those for 
the airplane. However, the horizontal -tail bending moments, perhaps fortuitously, are in 
good agreement over the test range of Cn- At a nominal M = 0.90 (fig. 12(f)), the pre- 
dicted wing buffet bending moments for the model are seen to scatter about the almost 
linear airplane values. Little can be said of the horizontal -tail data other than that both 
the model and airplane data exhibit considerable scatter. The c.g. buffet accelerations 
predicted from model data follow the trend of increasing severity with increasing Cjq, 
whereas the airplane measurements indicate little if any increase with Cjj. One might 
question the almost instantaneous increase in airplane buffet acceleration from the onset 
value of (presumably) a little over 0.3 percent of the design load factor to a value in excess 
of 0.7 percent of design. 

A = 72° .- At M M = 0.52 (fig. 12(g)), the airplane c.g. buffet accelerations are seen 
to remain constant with increasing C^, whereas the predicted values tend to increase 
initially with increasing to a maximum value slightly greater than the flight value 

and then to decrease slightly. The same trend is exhibited in the model wing and 
horizontal -tail bending moments. One may imagine that this trend is observable in the 
airplane wing buffet bending moments but no such trend is discernible in the airplane 
horizontal-tail data. No model data are presented for the right horizontal tail due to a 
failure of the strain gage. 
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At M m = 0.80 ^ = 0.81 to 0.76^ (fig. 12(h)), the model wing buffet bending 
moments again indicate increasing severity with increasing Cjg to an apparent maxi- 
mum valueandaregenerally larger, than Jhe flight buffet levels at the higher levels of 
C N . (Once again the airplane Mach number decreased substantially at the higher - normal- 
force values.) Both model and airplane horizontal-tail data show considerable scatter 
with the airplane bending moments generally greater than those predicted by the model. 
The trends of the model and airplane c.g. buffet acceleration are comparable although the 
model accelerations are higher than the airplane accelerations. Also, this is the first 
instance where the model c.g. acceleration variation with Cjj was not the same as the 
wing bending-moment variation. 

The limited model data at = 1.17 (fig. 12(i)) does not compare well with the 
airplane data. Model buffet onset occurred at a considerably lower value of than 

the airplane buffet onset, and lack of model tail power accounted for the small overlap of 
the model and airplane data. At comparable values of Cjq, the model and airplane wing 
buffet bending moments are in good agreement (no model data were obtained for the left 
wing due to strain-gage failure). Both the model and airplane exhibited very low wing 
buffet response (less than 1/4 percent of airplane design). The horizontal -tail bending 
moments predicted from model data are about the same magnitude as the airplane data, 
and the model c.g. buffet accelerations are relatively high compared with those on the 
airplane. 

It is felt that the model data at the higher sweep angles and higher Mach numbers 
could be improved by using a smooth skin type instead of the segmented type construction. 
The small discontinuities caused by the edges of the wing segments are more nearly par- 
allel to the airflow at the lower sweep angles but are presented more normal to the flow 
as the wing sweep increases. At the higher Mach numbers they provide sources for the 
generation of small shocks which can disturb the boundary layer and may cause premature 
separation. 

CONCLUDING REMARKS 

A wind-tunnel technique which employs complete dynamically scaled aeroelastic 
models for predicting airplane buffet loads has been evaluated. The buffet bending 
moments and accelerations measured on a flutter model of a fighter airplane are com- 
pared with those measured on the airplane in an independent flight buffet research study. 
Some pertinent observations follow. 

Although the model used in the study was a flutter model and was not specifically 
designed for buffet load studies, analysis of the model vibration characteristics and the 
modal character of the model and airplane structural response to buffet flow indicated 
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that the model adequately simulated the airplane structure. Even with this type of elabo- 
rate cable-mounted flutter model, practical considerations call for some compromises 
such as lack of simulation of pitch rates and Reynolds number. To improve boundary- 
layer simulation on the model, a smooth skin construction technique instead of the seg- 
mented approach is recommended. 

Potential improvements in technique were identified, such as providing supplemen- 
tary means of changing angle of attack so that horizontal -tail deflections could be made 
independent of model angle of attack. 

Power spectral density analyses of the model response to buffet flow showed that 
the indicated modal composition of the total measured response was dependent on the type 
and location of the measurement, wing sweep angle, Mach number, and, in some cases, 
the depth of penetration into the buffet region. 

The normal-force coefficients at which buffet onset occurred on the airplane were 
generally very close to those predicted by the model except for the 72° wing sweep con- 
figuration. At this sweep angle, buffeting began on the model at significantly lower 
normal-force coefficients than indicated for the airplane. 

The full-scale buffet bending moments on the wing and horizontal tails and the c.g. 
buffet accelerations predicted from the model data agreed well with airplane values at all 
Mach numbers at a wing sweep angle of 26°. At a wing sweep of 50°, the agreement was 
reasonably good at all Mach numbers tested for the wing bending moments, but the cor- 
relation of the model and airplane c.g. accelerations and horizontal-tail bending moments 
was not as good at the higher Mach numbers. At 72° sweep, both the airplane and model 
data exhibited a large degree of randomness at extremely low levels of buffet response 
which made evaluation of the correlation difficult. However, the airplane and model 
loads were of the same general magnitude. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., December 18, 1972. 
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TABLE I.- FLUTTER MODEL NOMINAL SCALE FACTORS 


Parameter 

Symbol 

Airplane-to-model ratio 

Reduced stiffness 

k r = (bo)/V) r 

*1.195 

Length 

b r 

8.0 

Mach number 

M r 

1.0 

Velocity 

V r 

*2.232 

Frequency 

U) r 

1/3 

Airstream density 

P r 

*1.457 

Dynamic pressure 

<lr 

*7.258 

Stiffness 

(EI) r , (GJ) r 

42 450 

Mass 

m r 

745.9 

Mass moment of inertia 

Ir 

47 736 

Mass density 

l-f-) 

*1.0 


\» H 



*For the buffet test conditions, these values vary slightly. (See section in text 
entitled "Scaling Considerations" and table III.) 
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TABLE II. - SAMPLE COMPARISONS OF STRUCTURAL DAMPING AND SCALED 
MODEL FREQUENCIES WITH AIRPLANE VALUES 


Symmetrical mode 

c s ; 

/C cr 

f, 

Hz 

Model 

Airplane 

Model 

Airplane 

A = 26 ° 

Wing first bending 

0.008 

0.009 

5.23 

5.15 

Wing second bending 

.012 

.010 

15.3 

15.3 

Wing torsion 

.010 

.020 

25.3 

25.4 

Fuselage vertical bending 

.008 

.009 

7.4 

7.6 

Horizontal -tail bending 

.014 

.010 

13.3 

13.3 

Horizontal -tail pitch 

.010 

— 

32.4 

32.7 

A = 50 ° 

Wing first bending 

0.007 

0.009 

4.9 

5.1 

Wing second bending 

.008 

.007 

15.1 

16.7 

Wing torsion 

— 

.018 

25.7 

26.4 

Fuselage vertical bending 

.005 

.009 

7.4 

7.5 

Horizontal -tail bending 

.016 

.020 

13.4 

13.3 

Horizontal -tail pitch 

.010 

— 

31.7 

32.7 

> 

ii 

CO 

o 

Wing first bending 

0.012 

0.010 

4.9 

5.1 

Wing second bending 

.007 

.022 

15.2 

16.2 

Wing torsion 

— 

.020 

27.7 

26.6 

Fuselage vertical bending 

— 

.009 

7.4 

7.5 

Horizontal -tail bending 

— 

.018 

13.0 

13.3 

Horizontal -tail pitch 

— 

— 

31.0 

32.7 
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A, deg 

M 


D 


26 

0.52 

1.582 

2.130 

7.18 


0.76 

1.945 

1.923 

7.20 


0.81 

1.962 

1.970 

7.62 

50 

0.52 

1.618 

2.121 

7.29 


0.80 

2.075 

1.952 

7.90 


0.90 

2.120 

1.928 

7.86 

72 

0.52 

1.618 

2.103 

7.17 


0.80 

2.185 

1.952 

8.34 


1.17 

1.704 

2.050 

7.22 


k 1/2 
K r 

k d 

Mode 

(a) 


/ \ - 1 / 2 
fei 

K a 

1.119 

0.859 

W1B 


1.065 

4380 



FVB 

.1676 

1.072 

— 



HTB 

.1670 

1.085 

4455 

1.117 

0.962 

W1B 

0.0878 

BB 

4400 



FVB 

.1676 

mm 




HTB 

.1670 

mm 

4470 

1.164 

0.995 

W1B 


0.998 

4535 



FVB 


.956 

— 



HTB 

.1795 

1.013 

4600 

1.121 

0.883 

W1B 

0.1282 

BB 

4410 



FVB 

.2272 

Mm 

— 



HTB 

.2463 

1.052 

4410 

1.168 

1.043 

W1B 

0.1017 

0.972 

4590 



FVB 

.1802 


---- 



HTB 

.1955 


4580 

1.172 

1.052 

W1B 

0.0885 

0.968 

4235 



FVB 

.1568 

.965 

— 



HTB 

.1700 

.990 

4330 

1.125 

0.876 

W1B 

0.1270 


4420 



FVB 

.2403 





HTB 

.2138 


4365 

1.169 

1.097 

W1B 


0.968 

4835 



FVB 

.1823 

.949 

— 



HTB 

.1622 

.954 

4765 

1.140 

0.900 

W1B 

0. 1403 


4455 



FVB 

.2652 

1 

— 



HTB 

.2360 


4400 



0.738 


0.768 


0.762 


0.732 


0.795 


0.737 






































































































16 Hz, Wing first bending 



34 


(b) Variation of sidewall fluctuating pressure power spectrum at pertinent discrete frequencies as a function 
of Mach number in Langley transonic dynamics tunnel. Test medium, Freon-12; empty tunnel. 

Figure 1.- Concluded. 



flying cable 



L-68-3209.1 

mounted in Langley transonic dynamics tunnel. 










































Figure 4.- Lift-balancing device for testing models at high lift coefficient. 
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Measurements: 

Left and right wing root bending moments 

Left and right horizontal-tail root bending moments 

Dynamic acceleration at c.g. 
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Figure 5.- Schematic diagram of data recording instrumentation. 





















Normalized bending-moment power spectral density, 
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Airplane, Op = 1.065 percent design 
Model, Op = 1.043 percent design 

(D Wing first bending 


Airplane frequency, Hertz 
(a) Right wing bending moment. 

Figure 8.- Sample comparisons of model and airplane response spectra 
A = 26°; M = 0.81; spectra normalized on rms level. 



Bending-moment power spectral density. 


100 xlO 


Airplane, Op = 0.99 percent design 

Model, Op = 1.08 percent design 

(1) Wing first bending 

(2) Fuselage vertical bending 

( 3 ) Horizontal-tail bending 

( 4 ) Vertical -tail torsion 

( 5 ) Horizontal-tail pitch 


Airplane frequency, Hertz 
(b) Left horizontal-tail bending moment. 
Figure 8.- Continued. 






Airplane, = 0.85 percent design 

Model, a^ = 1.50 percent design 

(5) Fuselage vertical bending 

(2) Wing second bending 

( 3 ) Wing torsion 

( 4 ) Horizontal-tail pitch 


Airplane frequency,, Hertz 
(c) Acceleration at center of gravity. 
Figure 8 .- Concluded. 




, percent design ®r» percent design . percent design 


3 . 0 - 





(a) Total response at A = 26°. M M = 0.52; M A = 0.53 to 0.49. 


Figure 9.- Model and airplane total response illustrating method of determining 
levels of buffet response and onset values. 
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™ j percent design ®j>» percent design a>p, percent design 
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Figure 10.- Comparison of buffet acceleration spectra at center of gravity and 

pilot station. A = 26°; M = 0.52. 
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O c.g. 

□ Pilot station 


M * 0.52 









, percent design a A , percent design a . , percent design 


c.g. acceleration 

O Model 
• Airplane 


Indicated airplane 
buffet onset — _ 


«p • 


Horizontal -tail root bending moment 
Model Airplane 
O Right • 
a Left r 


° * i 
J° ° 

o c * • 

% 4 


Wing root bending moment 
Model Airplane 
O Right • 


(a) A = 26°; M M = 0.52; M A = 0.53 to 0.49. 

Figure 12.- Comparison of airplane buffet response and buffet response predicted 
from model data normalized on airplane design loads. 
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(g) A = 72°; M m = 0.52; M a = 0.56 to 0.51. 
Figure 12.- Continued. 
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(i) A = 72°; M m = 1.17; M A = 1.23 to 1.16. 

Figure 12.- Concluded. 
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